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Abstract. The dual of a matrix ordered space has a natural ma- 
trix ordering that makes the dual space matrix ordered as well. 
The purpose of these notes is to give a condition that describes 
when the linear map taking a basis of M n to its dual basis is a 
complete order isomorphism. We exhibit "natural" orthonormal 
bases for M n such that this map is an order isomorphism, but not 
a complete order isomorphism. Included among such bases is the 
Pauli basis. Our results generalize the Choi matrix by giving con- 
ditions under which the role of the standard basis {-Ey} can be 
replaced by other bases. 




Given a vector space V there is no "natural" linear isomorphism 
between V and the dual space V d , but each time we fix a basis B = 
{vi : % G /} for V there is a dual basis B = {5i : i E 1} for V d satisfying 

6i(vj) = 

and this allows us to define a (basis dependent) linear isomorphism 
between V and V d . 

Definition 1. If B is a basis of M n , the linear map from M n to M d 
taking each member of B to the corresponding member of the dual 
basis is denoted by T>b, and is called the duality map. We let = 
Dg 1 : M d ^ M n denote the inverse of this map. 

Note that if / E M d , and B is a basis of M n , then T B (f) = J2 be ts f( b ) b - 
In particular, when V = M n (the space of n x n complex matrices), 
and we let S = {E it j : 1 < i, j < n} denote the standard matrix units, 
then the map T £ : M d — > M n satisfies 

n 

»(/) = J2 f( E ^r 
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Definition 2. If / G M%, there is a unique matrix D such that f(X) = 
tr(DX) for all X G M n , and we call this matrix the density matrix for 
/, with no requirement of positivity for / or D. 

Thus is just the map that identifies a functional / with the trans- 
pose of its density matrix: 

(1) f(X) = tr(r £ (/)*X) for all X G M n . 

This note is motivated by the following result of Paulsen- Todorov- 
Tomforde [31 Thm. 6.2], which we will see later is a restatement of the 
Choi-Jamiolkowski correspondences [21 H]. 

Theorem 3. The map T>£ : M n — > M% is a complete order isomor- 
phism between these matrix ordered spaces. 

In this paper, we will show that this theorem is very dependent on 
the choice of basis. In fact, we will show that there exist orthonormal 
bases for M n such that the inverse duality map does not even send 
positive functionals to positive matrices. Even more intriguing, we will 
show that there are "natural" orthonormal bases for M n such that the 
inverse duality map does send positive functionals to positive matrices, 
yet does not send completely positive maps to positive block matrices. 
These results can be interpreted as giving some new Choi-Jamiolkowski 
type results. 

We believe that such bases might be useful as entanglement wit- 



nesses. (We will comment further on this after Corollary 12). 

Before proceeding it will be necessary to establish some notation. 
Recall that when we say that a vector space V is matrix ordered we 
mean that for each natural number p, we have specified a cone C p in the 
vector space of p x p matrices over V, M P (V), which we identify as the 
positive elements in M p (V), and that these cones must satisfy certain 
natural axioms, such as if A G C p and B G C q , then A © B G C p+q . 
We also require that if X = (xij) is a p x q matrix of scalars and 



A = (vij) G C p , then 



X*AX = Xk,iVk,iXi,j) e C q 



k,l=l 



When there is no ambiguity we simply write C p = M P (V) + . See [HI 
dipt. 13] for more background on matrix ordered spaces. 

Matrix ordered spaces are the natural setting for studying completely 
positive maps. Indeed, given matrix ordered spaces V, W we say that 
a linear map $ : V — > W is completely positive provided that for each 
p, (vij) G M p (V) + implies that ($(^)) G M p (W) + . (We will denote 
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the map that takes {vij) to (^(fij)) by so that $ is completely 
positive iff every map $^ is a positive map.) We say that $ is a 
complete order isomorphism provided that $ is invertible and that $ 
and $ _1 are both completely positive. 

The most frequently encountered example of a matrix ordered space 
is L(H), the bounded linear operators on a Hilbert space H. We define 
the matrix ordering by identifying M P (L(H)) = L(H © ■ ■ • © H), the 
bounded linear operators on the direct sum of p copies of the Hilbert 
space, and declaring (Aij) e M P (L(H)) + exactly when it defines a 
positive operator on the Hilbert space H ffi • • ■ © H. 

When H = C n , the standard n-dimensional Hilbert space, we write 
ei, . . . , e n for the standard basis of C n , and write (■ , ■) for the inner 
product on C n . We identify L(C n ) with M n , the set of n x n matri- 
ces with entries in C. Note that identifying M p (M n ), the p x p block 
matrices with entries from M n , with M pn yields the same cone of posi- 
tive matrices as when we identify M p (M n ) with the linear maps on the 
direct sum of p copies of C n , L(C n © • • • © C n ). 

Finally, given a matrix ordered space V there is a natural way to 
define a matrix ordering on the dual space V d . To do this we declare 
that a matrix of functionals (f%j) G M p (V d ) belongs to M p (V d ) + if 
and only if the linear map $ : V — > M p given by $(t>) = (fij(v)) is 
completely positive. 

In this paper we will be concerned with examining various natural 
bases B for M n and determining whether or not the duality map Z)g is 
a complete order isomorphism. We will see that there exist bases for 
M n such that D$ is positive but not completely positive. 

Since our results rely on Theorem [3j we present a new proof here, 
that is somewhat simpler than the proof that appeared in [3]. 

Proof of Theorem [3| Rather than proving that T>£ is a complete 
order isomorphism, we prove, equivalently, that = V^ 1 : M d — > 
M n is a complete order isomorphism. We have already seen that T^- 
sends functionals to the transpose of their density matrices. Since a 
functional is positive if and only if its density matrix is positive, and 
the transpose map is an order isomorphism, we see that is an order 
isomorphism. 

Now let (fk,i) e M p (M d ) and consider the map $ : M n — >■ M p 
defined by $(X) = (f k ,i(X)) = £* l=1 f k ,i(X)E k j. We must show that 

$ is completely positive if and only if Ig ((/j^)) G M p (M n ) + . 
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We have 

n 

r i P) ((A,z)) = (^e(fk,i)Y kl=1 = ( ^2 fk,i(Eij)E id )) P kl=1 

p n n p 

fc,i=l 2ji=l »jjf=l A;,Z=1 

n 

= ^ '!'( /•.;.,) /•-•;..;• 

Since the map that takes A® B to B ® A extends to a ^isomorphism 
of M p £g> M n onto M n (g) M p , the last expression is positive iff the matrix 

n 

(2) C„ = Y E a ® 

is positive. But this last matrix is the Choi matrix and by Choi's 
theorem [2] the map $ is completely positive if and only if this block 
matrix is positive. Thus (fk,i) G M p (M n ) + if and only if (Ts(fk,i)) £ 
M p (M n ) + and we have shown that Tf is a complete order isomorphism. 
This completes the proof of Theorem 1. □ 

A map \I/ : M n — )■ M n is called a co-positive order isomorphism 
provided that its composition t o with the transpose map t on M n is 
a complete order isomorphism. 

Corollary 4. TTie linear map from M% to M n that takes a functional 
to its density matrix is a co-positive order isomorphism. 

Proof. As remarked in connection with equation (JTJ, the map that takes 
a functional to its density matrix is t o Tg. □ 

Now that we have a complete order isomorphism T>g between M n 
and M„, in order to determine whether or not other maps between M n 
and M% are complete order isomorphisms, it will be convenient to work 
with a map V B : M n -> M n instead of V B : M n -> M%. 

Definition 5. Let B be a basis of M n and 8 the standard basis of 
matrix units. Then we define T> B : M n — > M n by T> B = Tg o T>%. 

Note that since Tg is a complete order isomorphism, T> B will be 
a complete order isomorphism if and only if T> B is a complete order 
isomorphism. 
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Given a matrix A G M n we let A 1 denote its transpose. Recall that 
£ = {Eij} is an orthonormal basis for M n with respect to the Hilbert- 
Schmidt inner product, which we denote by (A,B) = tr(AB*), where 
tr : M n — > C denotes the unnormalized trace, tr(A) = Y^i+i °M- 

Using the orthonormal basis £ for M n , we represent elements in 
L(M n ) as n 2 x n 2 matrices. Given L G L(M n ) we write L T to denote the 
transpose of the matrix for L with respect to the basis £, to distinguish 
this transpose from the transpose on M n . 

Definition 6. Let B be a basis of M n and £ the standard basis of 
matrix units. A change of basis map is any linear map Cg in L(M n ) 
taking the set £ to the set B. By slight abuse of notation, we write 
Cq for the unique linear map in L(M n ) whose matrix in the standard 
basis £ is the transpose of the matrix of C#. We define Mg = CbCq G 
L(M n ). 

Since a linear map is uniquely determined by its values on a basis, 
we see that a change of basis map is uniquely determined up to re- 
orderings of the basis elements. Thus, in the setting of M n there will 
be (n 2 )! change of basis maps. However, the map Mg is independent 
of the choice of change of basis map or matrix. Indeed, fix one change 
of basis map Cg. Then every change of basis map has the form o P, 
where P G M n is a linear map which permutes the basis £. Since 
P sends the orthonormal basis £ to itself, its matrix in that basis is 
orthogonal. Thus M B is unchanged if we replace Cq by CqP. 

Fortunately, we will find that the results that we seek depend on 
Ms and are independent of the choice of change of basis map Cg. In 
particular, our conditions determining when the duality map T>b is an 
order isomorphism or a complete order isomorphism will be expressed 
in terms of the map Mg. 

Theorem 7. If B is a basis of M n , then the duality map is given by 
V B = V £ oM B \ 

Proof. Let B = {X u . . . , X n n}, and let {X u . . . , X„ 2 } c M% be the 
dual basis. We write £ = {E\, . . . , E n 2}. 

Define Y 6 = V^ 1 (X j ) = T e (Xj) G M n . We are going to show 

(3) C^Yj = Ej for all j. 
For each i, j, using ([T]) 

(E^C^Yj) = (CgE^Yj) = (X^Yj) = tr(X,F*) = tr(XY/) 

(4) =tr(X l (T E (X J )) t )=X J (X, i )=5 ir 
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It follows that CgYj = Ej. Taking conjugates gives 
Now 

CeCgD^Xj = CbCqYj = C B Ej = Xj = D^Xj, 
so by linearity M B D £ l = D B l . Thus D B = D £ M B l . 

□ 

Notation. If C G M n , then $ c : M n — > M n is the completely positive 
map defined by $ C {X) = CXC*. 

Recall that a map \I/ : M n — >■ M p is called completely co-positive if 
and only if t o ^ is completely positive, and a map \I/ : M n — > M n is 
called a co-positive order isomorphism provided that its composition 
t o with the transpose map t on M n is a complete order isomorphism. 
Here the order of composition with the transpose map doesn't matter 
as shown by the next result. 

Proposition 8. Let $ : M n — > M p then t o $ is completely positive if 
and only if $ o t is completely positive. 

Proof. By Choi's result [2J, t o $ is completely positive if and only if 
[^(EijY) is positive, while $ o t is completely positive if and only if 
(&(Ej t i)) is positive. But these np x np block matrices are transposes 
of each other. □ 

We start with the following description of order automorphisms of 
M n and their partition into completely positive and completely co- 
positive maps. It is a consequence of more general results of Kadison [5] 
relating isometries, Jordan isomorphisms, order isomorphisms, and *- 
isomorphisms of C*-algebras, specialized to M n viewed as a C*-algebra. 

Lemma 9. Let $ : M n — > M n be an order isomorphism. Then there 
exists an invertible C G M n such that either $ = $c* or $ = t o <& c . 
In the first case, $ is a complete order isomorphism, and in the second 
case $ is a co-positive order isomorphism. If n > 1, both cases cannot 
occur simultaneously. 

Proof. First assume $ is unital, i.e., that = /. For Hermitian 

matrices A, we have \\A\\ = sup {A G R | —XI < A < XI}, so a 
unital order isomorphism is an isometry on Hermitian elements of M n . 
It follows that $ is an isometry on all of M n (5J proof of Thm. 5]. 
Every unital isometry on M n is a Jordan isomorphism, i.e., preserves 
the Jordan product A o B = (1/2) (AB + BA), cf. Thm. 7]. Every 
Jordan isomorphism on M n is a *-isomorphism or *-anti-isomorphism 
[SI Cor. 11]. In the latter case, composing with the transpose map 
gives a *-isomorphism. It is well known that every ^isomorphism of 
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M n is conjugation by a unitary, see, for example, [U Thm. 4.27]. Thus 
there is a unitary U such that $ = $[/ or $ = t o 

Finally, let $ be an arbitrary order isomorphism. We will show 
$(/) is invertible. Observe that < A G M n is invertible iff A is an 
order unit, i.e., if for all B = B* G M n there exists A 6 R such that 
— XA < B < XA. An order isomorphism takes order units to order 
units, so $(1) is invertible. Let D = $(1) 1/2 , and define * = $ D -i o $. 
Then ^ is a unital order isomorphism, so by the first paragraph there 
exists a unitary U such that \I/ = $[/ or \I/ = toQjj. Then $ = $£>o\I> = 
$d ° &u = &du, or else $ = $ D oto$ l/ = to$ c . where C = D t U. 
Note that if both cases occur then \1/ is both a ^isomorphism and a 
*-anti-isomorphism, which is possible only if M n is abelian, and that 
holds only when n — 1. □ 

Note that Lemma [9] implies that a composition of two complete order 
isomorphisms, or two co-positive order isomorphisms, is a complete 
order isomorphism, and a composition of a complete order isomorphism 
and a co-positive order isomorphism (in either order) is a co-positive 
order isomorphism. 

Theorem 10. Let B be a basis of M n . Then T>g is an order isomor- 
phism iff there exists C 6 M n such that either (1) Mb = $c or (2) 
Mg —to <& c . In the first case, a complete order isomorphism, 

and in the second case it is a co-positive order isomorphism. 

Proof. By Theorem [7| we have = T>£ o M^ 1 . Since is a com- 
plete order isomorphism (Theorem [3]), then Dq is a complete order 
isomorphism (respectively, co-positive order isomorphism) if and only 
if Mg is a complete order isomorphism (respectively, co-positive order 
isomorphism) . Now the theorem follows from Lemma |9} □ 

Corollary 11. Let B = {Bj : 1 < j < n 2 } be a basis for M n . 

(1) M B = $c for some C G M n if and only if T B : M% -> M n is a 
complete order isomorphism, 

(2) Ml = t o $ c for some C G M n if and only iftoT B : M^ — )■ M n 
is a complete order isomorphism. 

Proof. We prove the second statement. By Theorem [lUJ Mg = toQ c iff 
Djg is a co-positive order isomorphism. This is equivalent to Tg = Dg 1 
being a co-positive order isomorphism, and hence to t o Tg being a 
complete order isomorphism. 

□ 

Both Choi and Jamiolkowski have defined useful correspondences 
that associate a matrix in M n ®M p with each linear map $ : M n — > M p . 



8 



V. I. PAULSEN AND FRED SHULTZ 



Choi's correspondence is $ i— > C$, where 

(5) C* = ®$(%). 

y 

As remarked in the proof of Theorem |3j positivity of the Choi matrix 
^ is equivalent to positivity of £\ . ® E 1 ^-, and it is this latter 

form that we generalize below. We now describe a related correspon- 
dence defined by Jamiolkowski [4J. If $ : M n -»■ M p , then J"($) is 
defined by the condition ( J"($), A* ® B) = ($(A), B) for all A G M n , 
-B G M p . This is equivalent to 

(6) J($) = 5^J5ft®*(£k). 

Regarding our current investigation, the Choi matrix C$ has the prop- 
erty that C$ > iff <J> is completely positive, j2]. The Jamiolkowski 
correspondence has the property that in (|6]), J7"($) is unchanged if the 
basis {-Ejj} is replaced by any orthonormal basis of M n . Our corre- 
spondence will be closer to Choi's. 

Corollary 12. Let B = {Bj : 1 < j < n 2 } be a basis for M n , and let 
\& : M n — » M p be a linear map. 

(1) If Mb = $c f or some C G M n , then $ zs completely positive if 
and only if Y,f =1 V(Bj) <g> G (M p g> M n ) + . 

(2) If Ms = to $ c /or some C G M n , i/ien ^ is completely positive 

if and only if *( B j) ® B j e ( M p ® M n) + ■ 

(3) If Mq = t o <& c /or some C G M n , t/ien ^ z's completely co- 

2 

positive if and only if YTj=i ^i B j) ® B j e {M p <8> M n ) + . 

Proof. To prove the first statement, for 1 < k, I < p define G 
by ^(X) = (fk t i(X)). Then by the definition of the order on 
M P (M£) discussed earlier, ^ is completely positive if and only if (fk,i) G 
M p (M%)+ which holds if and only if (T B (f h> i)) G M p (M n ) + . But as in 
the proof of Theorem 1, we have that 

n 2 

(TbC/jm)) = X>(*i) ® s i- 

i=l 

To prove the second statement, note that \I/ is completely positive 
if and only if (t o T B (f k i)) G M p (M n ) + and this matrix is seen to be 
equal to 

n 2 

i=i 
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For the third statement, replace \I/ by £ o v|> in the second statement. 

2 

This shows \l/ is completely co-positive iff Xy=i ^(BjY ® £ (M p ® 
M n ) + , and now applying the transpose map gives (3). □ 

We now point out that bases with the properties indicated in the 
Corollary are related to entanglement witnesses. Indeed suppose B = 



{Bi, . . . , B n i\ is a basis of M n for which Corollary 12 (1) holds. Taking 
$ = / we have £\ Bi <g> B { > 0. Let $ : M n ->■ M n be a map that is 
positive but not completely positive. Define 

B = Bi <g> Bi and B® = ^ $(-Bj) g) 

i i 

Then for any positive X, Y, since $ > 0, we have 

(£$, x ® y> = $(5,) ® 5,, x ® y) 

i 

(7) = (($ ® i)5 , x (8> y) = (B„, $*(X) ® y) > 0, 

and hence 5$ is > on all separable states. Since $ is not completely 
positive, then 5$ ^ 0, so there is a state A such that A) ^ 0. 
Such a state is then entangled, so B<$> is an entanglement witness. 

Examples 

Notation. If x,y G C n , then R x>y G M n is the rank one operator 
defined by R x>y z = (z,y)x. Observe that £y = R ei , ej - 

Proposition 13. Let (Ay) G M n , with all Ay nonzero, and let B be the 

basis {XijEij}. Then T> B is an order isomorphism if and only if the ma- 
trix (A? ) is positive semi-definite with rank one. In that case, there are 
scalars a\, . . . , a n such that Ay = a^aj, and if C = diag(«i, . . . , a n ), 
then T>q = $c ; and hence T>b is a complete order isomorphism. 

Proof. Note that Cg is diagonal for the standard basis of M n , so Cg = 
C B . Thus Mjs(Eij) = (C B Cl)(E l3 ) = A|.£y. 

Suppose that the map : M n — > M n is an order isomorphism. We 
first consider the case where V^ 1 = M B = $ c . Then $ c -(£'y) = A^-Ey, 
so 

The ranges of the two sides must coincide, so for each i there is a scalar 
ati such that Cti = a^i. Substituting into the displayed equation gives 
XfjEij = aiOjEij, so Afj = aiaj. Thus the matrix (A?-) has rank 
one and is positive. Conversely, if (Ay) is positive with rank one, 
then there are nonzero scalars ai,...,a n such that aia] = Ay. If 
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C = diag(«i, . . . ,a n ), then one readily verifies that Mg = $c. Then 
"Dg 1 is a complete order isomorphism, and hence so is X>g. 
Now we examine the possibility that CqCq = to$ c . Then 

Let D = C*\ so that A?£y = DE jt D*. Then for all i,j 



^HjRe^ej — RDej,D ei - 



This implies that Dej is a multiple of ej for all i, j, which is impossible. 

□ 

Example 14. If C G M n is invertible, then = $ C t, so &c®c = 
^cc*- Hence for the basis B = {$c(Eij)}, we have Mg = &cch so by 



Theorem 10 B has the property that the map from this basis to its 
dual basis is a complete order isomorphism. In particular, if {-FV,} is 
any system of matrix units for M n , there is a unitary V such that <&y 
satisfies Qy{Eij) = Fy for all i,j, and so the map from {Fij} to its 
dual basis is a complete order isomorphism. 

On the other hand, if U : M n — > M n is unitary with respect to the 
Hilbert Schmidt inner product and takes £ to a basis B, it need not 
be the case that the duality map V B is a complete order isomorphism, 



as can be seen from Proposition 13 with An = i and \j = 1 for 



(hj) (1)1)- Hence, not every orthonormal basis of M n has the 
property that the duality map is an order isomorphism. 

For our next application we study the Pauli spin matrices. 

Theorem 15. Let B = {o"o, o~i, on 03} be the Pauli spin matrices, i.e., 

1 0\ A) l\ A) -A A 



a ° =l v o ij' ai= ^i oj' a2= ^ oj' - 1 

T/ien i/je duality map T>q is a co-positive order isomorphism. Further- 
more, let ^> : M 2 n — > M p be a linear map. Then ^ is completely positive 
if and only if 

3 

n,...,j„=o 

z's a positive 2 n p x 2 n p matrix. 

Similarly, \I/ z's completely co-positive if and only if 

3 

^2 ® • • • ® O ® °ii ® ■ ■ • ® Vin 

n,...,i„=o 

z's a positive 2 n p x 2 n p matrix. 
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[Cb][C t b ] 



Proof. Let Cb be the linear map such that 

C-b(Eu) = cr , C B (E 12 ) = 0i, C B (E 2 i) = 02, C-b(E 22 ) = 03- 

Then the matrix for Mg = C&Cg in the standard basis of M 2 is 

(2 0\ 
2 
2 
\0 2/ 

which is twice the matrix of the transpose map t : M 2 — >• M 2 . Thus 
T>b in this case is a co-positive order isomorphism. 

Applying Corollary[l2j we see that a map \l/ : M 2 — > M p is completely 
positive if and only if 

3 

® 0j e (M p <g> M 2 )+. 

i=Q 

Using the explicit form of the Pauli matrices, we obtain that \1/ is 
completely positive if and only if 

"#(00) + #(03), ^(0i) + ^(0 2 ) 
^(01) - iV(a 2 ), ^(00) - ^(03) 

is positive in M 2 (M P ), which is identical to Choi's theorem. 

On M 2 n the tensored spin matrices B® n = {0^ <8> • • • <E> 0j n 
3} form an orthonormal basis which we will call the spin basis. The 
standard basis of matrix units of M 2 n consists of the tensor products 
of the matrix units of M 2 . The map Cg®n : M 2 n — y M 2 „ taking the 
standard basis of matrix units to this spin basis is then the tensor 
product of the maps on each factor M 2 , so (Cg®n)(CB®n) T will be the 
transpose map on M 2 n. Thus the map from the spin basis on M 2 « to 
its dual basis will also be a co-positive order isomorphism. 



< ij < 



Again applying Corollary 12 yields that a map, \1/ : M 2 n — >• M„ is 



completely positive if and only if 



<t0®4 <8>---<g>^. 



J2 

ii,...,i n =0 

is a positive 2 n p x 2™p matrix. 

Similarly, \1/ is completely co-positive if and only if 

3 

®---®o®°ii ®---®°^ 

ii,...,i n =0 

is a positive 2 n p x 2 n p matrix. 
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□ 

For our final application we study the map from the Weyl basis to 
its dual basis. We will compute the duality map for the Weyl basis, 
with the conclusion that this map is a complete order isomorphism for 
n = 2, but is not an order isomorphism for n > 2. Below n > 1 is a 
positive integer, and all indices are viewed as members of Z n . 

Definition 16. Let eo, . . . , e n -\ be the standard basis of C n , and B = 
{E ab | a, b G Z n } the corresponding matrix units. Let U, V G M n be 
defined by Vej = z^ej and Uej = e J+ i where z = exp(2iii/n) Then 
{^U a V b | a, b G Z n } is an orthonormal basis for M n which we call the 
Weyl basis W. 

The unitary matrices {U a V b \ a, b G Z n } are usually called the dis- 
crete Weyl matrices or the generalized Pauli matrices. 

Lemma 17. Define C w G L(M n ) by C w (E ab ) = j K U a V b . With re- 
spect to the standard basis of matrix units, we have the following matrix 
entries for Cyy and Cy^C w : 

[Cw]ab,cd = Z db 5b+ C ,a 0,nd [CyyC w } a b,cd = $b-dSa,c~2d- 

Proof. We have U a V b ej = z bj E j+a jej so 

U a V b = J2^E J+a ,y 



Thus 



ab,cd {Cw(Ecd), E a b) 

' (U c V d ,E ab ) 



71 



z dj E j+C j, E ab ) 



— ^2z dj 5 j+c , a 5 jtl 



n 

3 

(8) = i d Vc,a 

Now 

ab,cd abjk cd,jk- 

jk jk 

In the first factor [Cw] a &,jfc °f the last sum we use the expression Q for 
[Cw]a6, C (i with the substitutions c—tj and d k. In the second factor 
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[Cw]cd,jk we use ([8]) with the substitutions a — > c, b — > d, c — > j, d — > k. 
We get 

(C\\)C^y)ab,cd = ~ ^ ] zkb 3b+j,a zkd 3d+j,c- 



n 

jk 



The summands will be nonzero if and only i{j = a — b = c — d (mod 
n). Thus 



{CyvCw)ab,cd — —^a-b,c-d 



z k(b+d) 



n 

k 



The sum will be zero unless b + d = 0, in which case it has the value 
n. Thus 

(Cy\;Cw)ab,cd = Sa-b,c-d$b+dfl = 0~b-d$a,c-2d- 

□ 



Note that Lemma [T7| gives 
(9) (CwC^ v )(E cd ) = E c ^2d-di 

so in particular CwC^, acts as a permutation on the basis of matrix 
units. 

Corollary 18. For the Weyl basis W, the duality map T>w is a com- 
plete order isomorphism if n = 2, and is not an order isomorphism for 
n>2. 

Proof. If n = 2, then from ^ CwC^ is the identity map, and hence 
"Dyy is a complete order isomorphism. 

Now let n > 2. Suppose first (to reach a contradiction) that CyyCyy = 
$c for some invertible C £ L(M n ). Then by ([9]), 

E-d,~d = (CwCw)(Edd) = &cEdd = Rce d ,Ce d - 

Thus for all d there are scalars A^ of modulus one such that Cea = 
X d e- d . Then 

E c ~2d,-d = (Cy\;Cy V )(E cc i) = &c(E c d) = Rce c ,Ce d = ^c^dE-c-d- 

This implies c — 2d = — c, so 2c = 2d mod n for all c, d. This is 
impossible for n > 2. 

Now suppose CyvC^ = to $ c . Then again applying ([9]), 

E-d,-d = {CwCw){Ed,d) = ($c(Ed,d)y = (Rce d ,Ce d Y = (Rce d ,Ce d )- 



This implies that for all d, Ced is a multiple of e-d, and hence Ce d is 
a multiple of e~d = z~d- As above 

E c -2d,-d = {CwCw){Ecd) = (Rce c ,Ce d Y = RcT d fie~ c = A c ArfE__ c 

This implies c — 2d = —c for all c, d, which again is impossible for 
n > 2. □ 
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Remark 19. The Weyl basis for n = 2 is slightly different than the 
Pauli spin basis. Indeed, one has 

s/2 s/2\0 V' s/2 s/2 V °)' 



s/f s/2 \0 -i) ' V2 UV ' s/2 VI 
These are the Pauli spin matrices except for normalization and a miss- 
ing factor of i in the last. Again, applying Corollary [12] yields the usual 
Choi condition. 

Recall that for the Pauli spin matrices we found CgCj was the trans- 
pose map, so T>b in that case was co-positive. 

Remark 20. Note that for if we take tensors of the Weyl basis 
for M2, then we will get another basis for which the duality map is a 
complete order isomorphism. 

The Conjugate Linear Duality Map 

Using the fact that M n is a Hilbert space, we also have a canonical 
conjugate linear isomorphism between M n and its dual space. This 
map is unaffected by whether we make the inner product conjugate 
linear in the first or second variable, so we use the physics convention 
that inner products are conjugate linear in the first variable. Thus, the 
inner product on M n can be given by 

(A, B) = Tr(A*B) 

and the conjugate linear Hilbert space duality map is given by 

V d : M n -»■ M% where V d (A)(B) = Tr(A*B). 

The inverse of this map 

T d = V~ d l : M d n -> M n 

sends the linear functional Ja{B) = Tr(A*B) to the matrix A which is 
the adjoint of the density matrix. 

Proposition 21. The duality maps T> d and T d are conjugate linear 
complete order isomorphisms. 

Proof. By Corollary [4] the linear map that sends a functional to its 
density matrix is a co-positive order isomorphism. Hence, the map 
that sends a functional to the transpose of its density matrix is a com- 
plete order isomorphism. But a matrix (qj) is positive if and only if 
the matrix (cy) is positive. Thus, the duality map T d , which sends 
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a functional to the adjoint of its density matrix is a complete order 
isomorphism. Consequently, so is its inverse T>d- D 

The above result has a nice interpretation in terms of bases. Choi's 
characterization says that a map $ : M n — > M p is completely positive 
iff the matrix C$ defined in ^ is positive. As observed in Example 
14, in the definition of C$, the basis {E^} can't be replaced by an 



arbitrary orthonormal basis. The following result provides an alter- 
nate description of the Choi matrix that does have this independence 
property. Given a matrix B = (fry) we set B = (fojj). 

Proposition 22. Let {Bi}f^ be an orthonormal basis for M n . A com- 
plex linear map $ : M n — > M p is completely positive if and only if 

n 2 

(io) 

1=1 



is a positive np x np matrix. The matrix in (10) is independent of the 
choice of orthonormal basis, and equals the Choi matrix C$ . 

Proof. Let f A G M% be given by f A (B) = Tr(A*B), so that 

n 2 n 2 

F d (f A ) = A = J2(Bi,A)B l = YMBt)B h 
i=i i=i 
Thus, with respect to this basis 

n 2 

ra(/) = ]T7WB| for all f e M d n . 
i=i 

Let (f i:j ) e M% be the matrix defined by $(A) = (/# (A)) for A e M n . 
Recall that by definition, the matrix (/y) is positive iff $ is completely 
positive. Using the fact that is a complete order isomorphism, we 
have that $ is completely positive if and only if 

n 2 n 2 

(ii) (r,/(/,,)) = (MBijZi) = J2^b7)^b 1 

i=i i=i 

is a positive np x np matrix. Using that fact that a matrix is positive 
if and only if its complex conjugate matrix is positive, the proposition 
now follows by applying the ^isomorphism that takes A £g> B to B £g> A. 
Finally, since the matrix (fij) € M% is determined by $, the matrix 
Bi <g> <&(Bi) is independent of the choice of orthonormal basis {B{\. 



For the standard basis {Ejj} the matrix (10) is just the Choi matrix 



and hence the matrix in ( 10 ) equals the Choi matrix for all orthonormal 



bases {B{\. □ 
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Note that the matrix in (10) is the partial transpose of the matrix 



J($) defined by Jamiolkowski, cf. (J6J). 
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